Onsager reciprocity relations without microscopic reversibility 
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' In tins paper wo sin >\v that ( lusa^or Ylac'nlup time reversal 
£Nj properties of thermodynamic fluctuations and Onsager reci- 
procity relations for transport coefficients can hold also if the 
microscopic dynamics is not reversible. This result is based 
> on the explicit construction of a class of conservative models 
which can be analysed rigorously. 

\ Fundamental contributions to the theory of irreversible 
£^ processes were the derivation of the reciprocal relations 



C. Landim 

IMPA, Estrada Dona Castorina 110 
J. Botanico, 22460 Rio de Janeiro RJ, Brasil 
and LAMS de I'Universite de Rouen, Faculte de Sciences 
BP 118,F-76134 Mont-Saint- Aignan Cedex, France 



in 



for transport coefficients in states deviating only slightly 



from equilibrium and the calculation of the most prob- 
""£*»able trajectory creating a fluctuation near equilibrium. 
^The first result was obtained by Onsager in 1931 and 
Otjhe second one by Onsager and Machlup |2j in 1953. The 
ft, c alculation of the most probable trajectory relies on the 
reciprocal relations which in turn are a consequence of 
' («J microscopic reversibility. It turns out that the trajectory 
f^ih question is just the time reversal of the most probable 
trajectory describing relaxation to equilibrium of a fluc- 
is^tuation. The latter is a solution to the hydrodynamical 
equations. 

& ' These topics have received a certain amount of atten- 
tion in the physical litterature in the course of the last 
forty years. No rigorous results have been however es- 
tablished. More recently, this subjet has been taken up 
in various papers attempting more rigorous approaches : 
H , pi , ||, in the context of so called interacting particle 
systems and |^] , @] in a context of deterministic dynam- 
ical systems. 

In Q| we discussed the following question: is micro- 
scopic reversibility a necessary condition for the validity 
of the Onsager and Onsager-Machlup results? The an- 
swer to this question is far from obvious because in going 
from the microscopic to the macroscopic scale a lot of 
information is lost and irreversibilities at a small scale 
may be erased when taking macroscopic averages. 

In [Q we have exhibited a class of microscopic nonre- 
versible stochastic dynamics for which the time reversal 
rule of Onsager-Machlup is still valid even for fluctua- 
tions very far from equilibrium. This class of dynamics 



concerns dissipative one component systems with a hy- 
drodynamic equation of gradient type. Therefore there 
is no Onsager reciprocity relation to verify. 

In this paper we present a class of nonreversible multi 
component conservative models giving rise at the macro- 
scopic level to nonlinear purely diffusive equations in the 
terminology of ||. 

The equations are of the following form 



a 

dtp = ^2 du ' { D (p) ' du 'P} 

i=l 



(1) 



where p(u, t) = (pi(u, t), . . . , p n (u, t)) is a vector standing 
for the densities of different kinds of particles and D is 
in general a nonsymmetric n x n matrix. 

Associated to our models there is an entropy functional 
S(p) that is written as the integral of a density s(p) : 
S(p) = I s(p(u))du. 

The Onsager coefficients are defined in this context by 



L(p) = D{p) ■ R(p) 



(2) 



where the matrix R is determined by the entropy density 
s(p(u)) in the following way 



d 2 



ii.j 



dp l (u)dp J {u) 



(3) 



which is by definition a symmetric matrix. Onsager's 
reciprocity relations mean that L is a symmetric matrix, 
a property which holds for our models. 

In the physical literature one usually proves the 
Onsager-Machlup time reversal property from Onsager 
reciprocity relations. In our approach, we follow the op- 
posite order : we obtain a large deviation functional from 
which we prove the Onsager-Machlup time reversal prop- 
erty and compute the entropy. This in turn allows us to 
prove Onsager reciprocity relations. Our results are not 
restricted to the neighborhood of the equilibrium. 

For simplicity, we shall restrict ourselves to one di- 
mensional two component models but all analysis can be 
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carried out for any space dimension or for any number 
of components. As in [Q, we consider periodic boundary 
conditions. The systems considered in this paper differ 
from those of Q due to the conservative character of the 
dynamics. 

We consider an interacting particle system that de- 
scribes the evolution of two types of particles on a lattice. 
The stochastic dynamics can be informally described as 
follows. Fix a nonnegative function g : IN — > M + such 
that <?(0) = < g(i) for i > 1 and a finite range mean 
zero transition probability p(-) on "Z Q2 y yp{y) = and 
p(x) = for \x\ large enough). We shall assume the jump 
rate g to be Lipschitz and to diverge at the boundary : 
\g(k + 1) — g(k) | < la and lim^oo g{k) = oo. If there are 
ki, i — 1,2, particles of type i at a site x of Z5, at rate 
p(y)g{ki + k2){ki/ki + £2} one particle of type i jumps 
from site x to x + y. This happens independently at each 
site. 

To define precisely the process, we introduce some no- 
tation. We shall consider particles evolving on the dis- 
crete one dimensional torus with N points, denoted by 
%n (the integers modulo N). Sites of Zn are denoted by 
x, y and the configurations by the greek letter 77 = (771, 772) 
so that rji(x) stands for the total number of particles of 
type i at site x for the configuration 77. 

The generator Ln of this Markov process acts on func- 
tions / as 

N 2 2 

w = yE E p(y) T r +y f ( 4 ) 

j=l x,ye% N 

where the addition in 7/.s means addition modulo N; the 
operators T^' y are defined by 

(Tf ,J 7) = rx(v)vj{x)[f(Vi' v ,'ite) - fivi,m)} 

with r x (rj) = g(t]i(x) + ^{x)) / {qx{x) + 772(2;)} and ( x > y 
is the configuration obtained from £ letting one particle 
jump from x to y : 

!((z) if z^x,y 
C(z) - 1 if z = x (5) 
({z) + l if z = y. 

The operators T%' v are defined in a similar way. 

This process has two conserved quantities : the to- 
tal number of 771 -particles and the total number of 772- 
particles. It is therefore expected that for each fixed 
density pi > there should exist an equlibrium state 
with global density of 77^ -particles equal to pi. 

To describe these equilibrium probability measures, for 
each ipi,f2 > 0, consider the product probability mea- 
sure v N n m on N^ N x ]N Zn defined by 

^, V2 {(^6; v(x) = ki,£{x) = fc 2 } 

Z{ Vuv > 2 ) g(ki+k 2 )\ kilk 3 \ 



for fci > and &2 > 0. In this formula Z(ipi,ip2) is a 
normalizing constant and g(k)l stands for g(l) ■ ■ ■ g(k). 
A simple computation shows that these measures are in- 
variant for the Markov process with generator Ln and 
reversible if and only if the transition probability p(-) is 
an even function. 

Define p l : M + x H+ -> R+ by p l ((p 1 ,(p 2 ) = 
E,,n [77^ (0)1 and set p* = p\ + P2- One can check that 
p* is a smooth strictly increasing function of ip\ + ip% . De- 
note by a = a(p*) the inverse of p* : a(-) — (p*)^ 1 (■). We 
have that (pi = (pi/pi + p2)a{pi + P2)- To keep notation 
simple we shall denote by b(p*) the function a(p*)/p*. In 
conclusion, for each fixed density {p\,p2) we obtained an 
invariant state with total density of Tii-partciles equal to 
Pi. We shall from now on fix a density p = (pi, p2)- 

Let us consider now the unit interval TT — [0, 1) with 
periodic boundary condition and functions 7$ : TT — > _ZR + , 
i = 1,2 with global density pi : f T "/i(u)du = pi. The 
main object of our study is the empirical density p N (t) = 
«W, M fW) : ' 

V?(t,x) = ± £ Vi(t,y)S(x-^). (7) 

If we denote by Q^ i72 the distribution law of the tragec- 
tories p N (t) when the initial measure is concentrated 
on a configuration pair (77^,77^) such that p N (0) — > 
(ji(u)du,"f2(u)du) as TV | °°j it is possible to show that 
Q71 72 converges weakly as N goes to infinity to the mea- 
sure concentrated on the path p(u, t) that is the unique 
solution of 

dtp = {a 2 /2)d u {D(p) ■ d u p] 
p(O,-)=7(0- 

In this formula a 2 = J2 y V 2 p{y) and D = D(p) is the 
nonsymmetric diffusion matrix given by 

D(p)=b(p)I + b'(p)J(p) , (9) 

where / is the identity and J{p) is the matrix with entries 
Ji,j(p) = Pi- 

The above result is a law of large numbers that shows 
that the empirical density in the limit of large N be- 
haves deterministically according to equation (Q). We 
can now ask what is the probability that our system fol- 
lows a trajectory different from the solution of (||) when 
N is large but not infinite. This probability is exponen- 
tially small in N and can be estimated using the methods 
of the theory of large deviations introduced for the sys- 
tems of interest in || and |l0| . The main idea consists 
in introducing a modified system for which the trajec- 
tory of interest (fluctuation) is typical being a solution 
of the corresponding hydrodynamic equation, and then 
comparing the two evolutions. For this purpose, for each 
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pair of smooth functions Hi = Hi(u,t), we consider the 
Markov process defined by the generator 

N 2 2 

L Hj = ^y e H J ((y +X )/N^-H J ( x /N,t) p{y)T ,,,+v f 

2 — ' 

3=1 

with p(-) and Tj v as previously defined. The function 
H can be interpreted as an external field. 

The deterministic equation satisfied by the empirical 
density is now 

f d t p = (a 2 /2)d u [D(p) -d u p) ~a 2 d u [b(p)A(p,H)} 
I P(0, -)=7(0 

(10) 

where A(p, H) is the vector with components Ai — 
pid u Hi. 

Given a function p{u 1 1) twice diffcrentiable with re- 
spect to u and once with respect to t and such that 
Jrjp Pi(u,t)du = pi this equation determines uniquely the 
field H = (H 1 ,H 2 ). 

The probability that the original system follows a tra- 
jectory different from a solution of (|^) can now be ex- 
pressed in terms of the field H . We introduce the large 
deviation functional 



i=1 Jo Jtt 



dub{p)p l (d u H l ) 2 . (11) 



Let Q be a set of trajectories in the interval [0, to]. The 
large fluctuation estimate asserts that 



'71.72 



where 



Io,t (G) = inf J 0jto (p) 
peg 



(12) 



(13) 



The sign ~ has to be interpreted as asymptotic equality 
of the logarithms. 

From the equations (|l2]), dl3|), one sees that to find 
the most probable trajectory that connects the equilib- 
rium p to a certain state "f(u) one has to find the p(u, t) 
that minimizes /-oc.o(/°) m the set Q 7 of all trajectories 
satisfying the boundary conditions 



lim p(u,t) = p , p(u,0) = 7(u) 

t — > — OO 



(14) 



It is now possible to prove, following the same approach 
of [f|, that the unique solution of our variational problem 
is the function p*(u,t) defined by 



p*(u,t) = p(u, -t) 



(15) 



where p(u, t) is the solution of the hydrodynamic equa- 
tion which relaxes to equilibrium with initial state 7. 



p* (u, t) is therefore a solution of the hydrodynamic equa- 
tion with inverted drift 



d t p=-(a 2 /2)d u {D(p)-d u p} . 



(16) 



Equation (|LJ) is the Onsager-Machlup time-reversal re- 
lation. 

Denote by 5(7) the functional defined by 



S(7) = inf I-oofiiP) 
peG-, 



(17) 



which, by the Boltzmann-Einstein relationship, has to be 
identified with the entropy of the system. By inserting 



|15j) in ( 17) we obtain an explicit formula for the entropy : 



s{^{u))du 



IT 



where 



2 

3=1 

and E{p) = J" log p'dp', F{p) = J p log b(p')dp'. It is 
possible to check that S(p(-,t)) decreases in time if p(-, t) 
is a solution of the hydrodynamic equation (|^) . 

Of course, the entropy could also be calculated from 
the equilibrium measure ^ and it is easy to see that the 
two expressions coincide up to an additive constant. 

This explicit expression for the entropy 5(0 permits to 
check Onsager's relations in our model. A simple compu- 
tation shows that the matrix R defined by equation (|^) 
is such that 



= 5.. 



1 



b'{l{u)) 

7i(w) b(j(u)) 



(18) 



where &y stands for the delta of Kronecker and 5(7) is 
the entropy density. The product L = DR can now be 
computed using the explicit formula for D given in (|^) 
and shown to be a symmetric matrix. 

In all the above calculations we never used the sym- 
metry properties of the transition probability p(-) so that 
they are valid both for reversible and irreversible dynam- 
ics. 

This provides conclusive evidence that macroscopic re- 
versibility, in the sense of validity of the above results, 
does not require microscopic reversibility. 
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